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Abstract. We investigate the horizontal distribution of zeros of the derivative of the Ric- 
mann zeta function and compare this to the radial distribution of zeros of the derivative of 
the characteristic polynomial of a random unitary matrix. Both cases show a surprising bi- 
modal distribution which has yet to be explained. We show by example that the bimodality 
is a general phenomenon. For the unitary matrix case we prove a conjecture of Mezzadri 
concerning the leading order behavior, and we show that the same follows from the random 
matrix conjectures for the zeros of the zeta function. 

<N 

i — i I. Introduction 

H 

The zeros of the derivative C'( s ) °f the Riemann zeta-function are intimately connected 
with the behavior of the zeros of ((s) itself. Indeed, a theorem by Speiser [TH] states that 
the Riemann Hypothesis (RH) is equivalent to Ci s ) having no zeros to the left of the crit- 
ical line. Thus, understanding of the properties the zeros of C'( s ) can provide important 
tools and insight into the study of RH. After Speiser's article this idea was explored by 
Berndt [lj and Spira [20], but not much progress was achieved until the work of Levinson 
and Montgomery [llj, who proved a quantitative refinement of Speiser's theorem. They 
showed that ((s) and C'( s ) have essentially the same number of zeros to the left of the crit- 
ical line a = Re(s) = |, and proved that as T — > oo, where T is the height on the critical 
line, a positive proportion of the zeros of C'( s ) ue m the region 
CN 1 log log T 

CT<- + (l + e )-^^, 6>0. (1.1) 

Consider the group of unitary matrices U(iV) with probably distribution given by Haar 
measure, which is the unique measure invariant under the left and right action of U(iV) on 
^ itself. Such a probability space is often known as the Circular Unitary Ensemble (CUE). Let 

A(z) be the characteristic polynomial of a matrix in the CUE. In recent years evidence has 
been accumulated suggesting that, in the limit as T — > oo, the local statistical properties 
of £(s) can be modeled by the characteristic polynomials of matrices in the CUE where 
N ~ log(T/27r). The connection between the Riemann zeta-function and characteristic 
polynomials is extensive; examples include the distribution of the zeros of C( s )> its value 
distribution and its moments. (For a series of review articles on the subject see P2] and 
references therein.) 

Assuming that random matrix theory (RMT) provides an accurate description of C( s )? 
the horizontal distribution of the zeros of Ci s ) i n proximity of the critical line should be the 
same as the radial distribution of the roots of A'(z) close to the unit circle. This idea was 
first developed by Mezzadri [13], who determined the distribution of the zeros of A'(z) that 
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are very far from the unit circle and conjectured the leading order term of the distribution 
very close to the unit circle. In this paper we prove his conjecture. We also perform an 
analogous calculation for the Riemann zeta-function and conjecturally find that the result 
agrees with the RMT model. In addition, we do numerical computations in both cases 
and find a surprising feature in the distribution of zeros of the derivative, namely that the 
probability distribution is bimodal. 

2. The zeros of C'( s )- 

We have mentioned that the main motivation for studying the zeros of the derivative of 
the Riemann zeta-function is its connection with RH. Indeed, Levinson and Montgomery's 
result is the basis for Levinson's method [10J, which Conrey j2] used to prove that at least 
40% of the zeros of the zeta-function are on the line o = |. 

Levinson's method involves estimating a weighted average of the zeros of C'( s ) to the left 
of 2 + a,/ logT for some fixed a > 0. Thus, zeros of ('(s) in the region | < a < \ + a/ logT 
are an inherent loss in Levinson's method. It would be useful to understand the magnitude 
of this loss. Alternatively, if we could find a lower bound for the number of zeros of £' in 
this region we could improve the estimate for the number of zeros on the critical line. 
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Figure 2.1. Zeros of CO) (dot), ('(s) (triangle), (square), and (C7C)'( S ) 
(star), with imaginary parts in the range 1015 < T < 1040. 



Figure 2.1 gives a representative example of the location of zeros of ('(s) and the relation- 
ship to the zeros of ((s) and various other derivatives of the zeta-function. It illustrates that 
zeros of C'( s ) close to the critical-line correspond to closely spaced zeros of ({s). We make 
this statement precise in Section [6j Also, a zero of ('(s) seems to be "missing" when zeros 
of ( are particularly far apart or when there are two successive large gaps. Indeed, there 
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can't be a zero of C'( s ) between every pair of zeros of ((s) because the density of zeros of 
C(s) is 2^1og(T/27r) while the density of zeros of ('(s) is =- log(T/47r). So on average there 
is a "missing" zero of Ci s ) m each T interval of width 27r/ log 2 w 9.06. 

Conrey and Ghosh [I] and subsequently Guo [7j improved Levinson and Montgomery's 



result (1.1) and showed that a positive proportion of the zeros of ('(s) are much closer to 
the line <J =\- Indeed, for any fixed a > 0, the region 

a- l ->^- (2.1) 
2 - logT v ; 

contains a positive proportion of the zeros. Soundararajan [IB] made further progress and 
introduced the functions 

m~(a) := liminf — - — - 1, (2.2a) 



0<7'<T 



m 



T^oo JVi(T) 



T->oo 



^ —2 log T 

0<Y<T 



where N\(T) is the number of zeros /3' + ij' of C'( s ) with < 7' < T. Soundararajan proved 
that m~(a) > for a > 2.6, and conjectured that 

m(a) = m _ (a) = m + (a). (2-3) 

He also conjectured that m(a) is continuous, that m(a) > for all a > 0, and that m(a) — > 1 
as a 4 00. Zhang [23], Feng [6], Garaev and Yildirim [8], and Ki [9] proved refinements 
of Soundararajan's results. In particular Feng [6] showed that m~(a) > for all a > 
unconditionally of RH but assuming a conjecture on the frequency of small gaps between 
consecutive critical zeros of ((s). 



3. Zeros of derivatives of polynomials and statement of results 

Suppose f(z) is a polynomial with all zeros on the unit circle. (Eventually, / will be a 
random polynomial obtained as the characteristic polynomial of a random unitary matrix.) 
The Gauss-Lucas theorem assures that all the roots of f'{z) lie on or inside the unit circle 
(zeros of /' on the circle occur only if / has multiple zeros). If f(z) has two zeros which 
are very close together, then f'{z) will have a zero close by. (This is a consequence of the 
continuous dependence of the zeros of /' on those of /. This dependence is actually piecewise 
analytic as will be described below.) The specific location of the nearby zero of f'{z) will 
depend primarily on how close those two zeros of f(z) are, and on the general position of 
the remaining zeros of f(z). Thus, to leading order (with respect to the size of the gap), 
the distribution of zeros of f'{z) near \z\ = 1 should largely depend on the distribution of 
small gaps between zeros of f(z), that is, on the tail of the nearest-neighbor spacing of zeros 
of f(z). We will make this idea precise and treat in detail the case where f(z) = A(z) is the 
characteristic polynomial of a CUE matrix (a matrix chosen from the unitary group XJ(N), 
uniformly with respect to Haar measure). 

Let z' be a root of A'(z) and define the random variable 



S := N(l - \z'\). 



(3.1) 
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Denote by Q(s;N) the probability density function (p.d.f.) of S. Mezzadri [13] showed that 
the limit 



Q(s) 



exists, and proved that 



Q(s;N) 



lim Q(s;N) 

Af->oo 



N — y oo, s — > oo, 



with s = o (N). He also conjectured that 

Q(s) 



,1/2 



3tt 



s -+ 0. 



(3.2) 
(3.3) 

(3.4) 



Formula (3.3) can be interpreted as the RMT counterpart of the Levinson-Montgomery 
bound ( 1.1 ) for the roots of ('{s). The RMT model of the Riemann-zeta function is based on 
the observation that the local correlations of the non-trivial zeros of ((s) coincide with those 
of the eigenvalues of matrices in the CUE. In order to make this correspondence quantitative, 
the densities of the eigenvalues and of the zeros of £(s) must be made (asymptotically) equal, 
i.e., 

N 
2^ 



1 T 
— log — . 



(3.5) 



It follows from (3.3) that the expected value of S does not exist — its p.d.f. does not decay 
sufficiently rapidly. On the other hand, using (3.3), the average of the values of S not 
exceeding N is 

nN ds 

s ■ — = logiV, iV -> oo. 



(3.6) 



Recalling the relation (|3.1|) between S and \z'\, we conclude that a positive proportion of the 

e from 

(1 + 6) 



roots of A'(z) must lie within a distance from the unit circle bounded from above by 

logiV 



N 



(3.7) 



from the unit circle. Because of (3.5), formula (3.7) corresponds to Levinson and Mont- 
gomery's result (1.1). 

Now consider the roots e ltl , . . . , e ltN (with — 7r < i; < tt) of the characteristic polyno- 
mial A(z) of a random unitary matrix distributed with Haar measure. It is convenient for 
our purposes to define 

xj = ^, j = l,...,N, (3.8) 

so that, on average, the distance between two consecutive x^s is one. The joint probability 
density function (j. p.d.f.) of the eigenvalues is given in terms of xi, . . . , xn as 

1 TT |2 



P2(Xi, 



.X N 



N N N\ 



where we have used the notation e^{x) 
1, . . . , N, if necessary, we assume that 



\ | |e7v(x fc ) - e N (xj)\ z , (3.9) 

l<j<k<N 

= ex.p(2nix / N) . Relabeling the indexes j 



X\ < . . . < xn < x\ + N. 



(3.10) 



ROOTS OF THE DERIVATIVE OF THE ZETA FUNCTION 5 

We also extend the sequence {xj} to be periodic by setting rcjv+i = xi, xn+2 = %2, etc. Fix 
integers n and j with < n < N — 2. Let us denote by pi{n\ s) the probability density 
function of Xj +n+ i — Xj. Since the j.p.d.f. (3.9) is invariant under translations, which means 



P2^xi + a, . . . , xn + ac) = P2(%i, • • • , xn) for all ael, (3.11) 

it follows that p2{n; s) does not depend on j. If n = 0, pz(s) : = P2(0; s) is known as the 
spacing distribution. It has an asymptotic expansion in powers of s: 

»w = i o - — i - (| - ^ + IB L J ) ,v 

1 +jAn'T^i)^^ + 0(s 7 ). (3.12) 



1 


1 


3 ~ 


3N 2 
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,315 



135iV 2 45iV 4 189iV 6 ^ 
This expansion follows from the pair correlation function for U(N) (see [15] 

sin 2 (7r?/) 



(3.13) 



N 2 sin 2 (ny/N)' 

and the fact that the pair correlation function and the nearest neighbor spacing for U(N) 
agree to order 6. 

To describe our main result, suppose that a root of A(z) is degenerate (which means 
j + i = Xj) so that z' = exp (2niXj/N) will also be a root of A'(z). Simple considerations 



,r 



of continuity show that, if Xj and slightly moved apart (say, while keeping the 

remaining roots of / fixed), then z' will also move, but will still be close to the midpoint 



of the segment joining exp(27uxj/iV) to exp(27rix-, + i/iV). In Proposition 8.1 we make this 
precise, showing that z' stays close to that midpoint provided Xj+i — Xj < l/n. Henceforth 
we assume that the rescaled distance 

6 := Xj +1 — Xj (3-14) 

is small. 

By the translation invariance of the j.p.d.f. of the Xj's, we may assume without loss of 
generality that 

6 6 
Xj + i = - and Xj = —-. (3.15) 

Define 

5:=N(l-z'). (3.16) 

In section |5] we shall show that 

5 = b 1 7T 2 6 2 + b 2 7c A 6 4 + 0{6 6 ), as 6^0, (3.17) 
where b\ and 62 are explicit functions of the zeros Xk for k ^ j,j + 1. 



By combining the distribution of 6 given in (3.12) with information we will determine 
about bi and 6 2 in Section [5j we will prove 

Theorem 3.1. Let A(z) be the characteristic polynomial of a random matrix in U(iV) dis- 
tributed with respect to Haar measure. The distribution of 5 = N(l — \z'\) arising from closely 
spaced zeros of A(z) is given by 

4 1/2 82 3/2 , n / 5/2 



— s ' s 

3tt 457r 

as N — > 00. 



+ 0(s 5/2 ), (3.18) 
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Note that Theorem 3.1 refers to the small values of 5 that arise from closely spaced zeros 



of the polynomial. The theorem does not account for all small values of S. That distinction 



is often missed, because examples such as shown in Figure |2.1| give the mistaken impression 
that zeros of the derivative very close to the unit circle can only arise from closely spaced 
zeros of the polynomial. Farmer and Ki [S] give examples of families of polynomials whose 
(rescaled) zeros are bounded away from each other, but for which the density function of 5 
vanishes like C- s as s — > 0. They also argue that any larger density of zeros of z' near the unit 
circle must arise from closely spaced zeros of the polynomial. Therefore we have the following 



corollary of Theorem 3.1, which proves Mezzadri's conjecture about the distribution of \z' 



Corollary 3.2. Let A(z) be the characteristic polynomial of a random matrix in U(iV) 
distributed with respect to Haar measure, and let Q(s; N) be the p.d.f. of S = N(l — \z'\) for 
z' a root ofA'(z). Then 

Q(s) = lim Q(s; N) = —s 1 ' 2 + 0(s). (3.19) 

AT-S.OO 37T 

Note that it remains an unsolved problem to show that Q(s) is a proper probability 
distribution. That is, to show J °° Q(s)ds = 1. This is the random matrix analogue of 
Soundararajan's conjecture m(a) — > 1 as a — > oo. 

4. Comparison with data 
We compare our formulas with numerical data. 

We generated Haar-random matrices in U(iV) using the simple algorithm described in 



Figure 4.1 shows the empirical distribution of the rescaled zeros of A' for various size ma- 



trices. Figure [42 shows the empirical cumulative distribution function Ip(x) = J* Q(s,40)ds 
for U(40) and a comparison with the tail of the empirical cumulative distribution function 
with our results, showing good agreement. 

We would like to know the underlying cause of the curious "second bump" in the dis- 
tribution of zeros of derivatives. This seems to be a completely general phenomenon. In 



Figure 4.3 we show the analogous distributions for characteristic polynomials of matrices 
from COE(40) and for degree-40 polynomials whose roots are independently and uniformly 
distributed on the unit circle. Both cases show the "second bump", although not quite at 



the same location. In Figure 6.1 we find a similar shape for the distribution of zeros of 



5. Proof of theorem 13.11 

Suppose f(z) is a degree- N polynomial having all zeros on the unit circle, for which two 
zeros Zi, z 2 are very close together. Then the derivative f'(z) will have a zero close to the 
midpoint {z\ + z 2 )/2. This follows because, if z\ = z 2 (that is, if f(z) has a multiple root at 
zi), then Z\ is also a root of the derivative f'(z), and the roots of /' are continuous functions 
of the roots of /. By a rotation we can assume that 

f( z ) = F(z)(z-e- ie / 2 )(z-e ie / 2 ), (5.1) 

where F(z) does not have any zeros e lt with —0/2 < t < 0/2. The root of /' near 1 is the 
root near 1 of 

jW = + T^hsr- + T^m' (5 ' 2) 

and we denote that root by z' = 1 — A. 
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Figure 4.1. Distribution of the (normalized) distance from the unit circle of 
zeros of A'(z) for A the characteristic polynomial of a random matrix from U(./V). 
Top row: N = 15, 40. Bottom row: iV = 100 and all three plots together. 
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Figure 4.2. Experimental cumulative distribution function of 5 for random ma- 
trices from U(40), and a comparison of the tail of the distribution with the main 
term in Theorem 13. 21 



We are concerned with the case when / is the characteristic polynomial of a random 
CUE matrix, and we want to understand the distribution of the zeros of /'. Since the 
CUE measure (i. e., normalized Haar measure on U(iV)) is invariant under rotation, the 
distribution depends only on the absolute value of the roots of /'. Those roots will accumulate 
near the unit circle as N grows, so we must rescale them suitably in order to get a meaningful 
result in the limit N — > oo. We let 

= 2n6/N 

A = SIN. (5 ' 3) 
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Figure 4.3. Distribution of the (normalized) absolute value of zeros of A'(z) for 
A the characteristic polynomial of a random matrix from COE(40) (left) and A 
a degree-40 polynomial with independent uniformly distributed (Poisson) roots on 
the unit circle (right). For the Poisson case, the plot is truncated to suppress the 
large contribution from the small values. 

Note that this rescaling gives (9) — 1. The rescaling of S is more subtle, and it will be found 



in equation (5.12) that this is the correct rescaling. Note that, while 9 is a real number, 5 is 



usually complex (but typically with small imaginary part, as we shall see). 

We will determine the leading order behavior of the roots which are close to the unit circle, 
so in the above notation we are interested in 



iV(l - \z'\) = N(l - |1 - A|) 



Re (5) - 



2N 



+ 



N 2 



(5.4) 



5.1. Expansion for the roots. Exploiting the symmetry of the j.p.d.f. (3.9) under arbitrary 
relabellings of the variables Xj, as well as its translation invariance (cf., Section 3.2), we will 



assume (without loss of generality) that F(z) from Equation (5.1 ) is given in the form (recall 
that eiy(x) := exp(27ra;/AQ) 

N-2 



F{z) = Yl (z - e N {x n )) 



x n e - 



n=l 



N 



u 



9 N 
2' 2~ 



(5.5) 



(We are not excluding the possibility that F has roots at ejv(±#/2).) Therefore, 

r-,; N-2 



^ ^ Z - e N{XnY 



Now let z = z' be a root of /' (say, that root which is closest to z — 1). The Implicit Function 
Theorem shows that z' is an analytic function of 9 (at least for 9 sufficiently small). Indeed, 
fixing F and regarding /' as a function f'{6; z) of both 9 and z, we have /'(0; 1) = 0; by the 
assumption that F has no root at z — 1 all that remains to observe is that 

0^2F(1) = jV(0;l) = r(0;l). (5.6) 



if is sufficiently large, by Proposition |8.1| z' is defined uniquely and analytically as a 
function of 9 in the domain 



\9\ < min{i, . . . , |xtv_ 2 |}. 



(5.7) 
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We write z' = 1 — 5/N and wish to expand (1 — 5/N — e7v(x n )) _1 as a Taylor series in 5. 
This will be justified when \S\ < N\l — ew(x n )\. Hence, we have 

F V 



j=0 



(5. 



for 5 sufficiently small, where 



JV-2 



F' 



U) 



(i) 



n=l 



e N (x n )) j+1 ' 



(5.9) 



We will see later that the prefactor of N is the right choice to make the coefficients A, 
approximately bounded. Note that, for \9\ < 1 (equivalently, for |0| < 2n/N), any 5 such 
that \5\ < \9\ = N\Q\/(2n) makes the expansion in ( |5.8 ) valid (independently of the exact 
location of the zeros of F), since 



U|<^<JNTI1 



2tt 



1 < n < N -2. 



e N (e/2)\<N\l-e N (x n )\, 
We also wish to expand the other terms in ( |5.2[ ) as a series in 5. We have 
1 1 2-f-2cos(f) 



(5.10) 



-i0/2 



+ 



oiB/2 



■25 

N 



+ ^-2(l-#) 



cos 



7t0_ 

N 



-25 + Ti 2 6 2 N- 1 - ±9 4 n 4 N~ 3 + 0(9 6 N 



-5> 



5 2 + tt 2 9 2 - jrWiV" 1 - ±tt 4 9 a N- 2 + -W 4 50 4 iV- 3 + 0(9 6 N~ 4 ) ' 



(5.11) 



Combining this with equations (5.2 ) and (5.8 ), putting all terms over a common denominator, 
and using the fact that f'{z') = 0, we have 



12 



7T 4 6 4 N- 2 + ±7r 4 56 4 N- 3 + 0(9 6 N~ 4 )) 
25 + vr^iV- 1 



= Aj5 j {5 2 + tt 2 9 2 - n 2 56 2 N- 1 

j=0 

-^TT 4 e 4 N- 3 + 0(6 6 N- 5 ). (5.12) 

Note that a global factor of N canceled to give the above equation, suggesting that we have 
chosen the correct scaling for 5. 

Equation (5.12) is simply a more explicit and manageable form of the equation f'(9; z') = 
defining z' implicitly as a function of 9. Noting that /'(0; 1) = 0, together with the functional 
equation f'(9; z) = f'(—6; z), it follows that 5 = 5(9) has an expansion in powers of 9 2 , with 
no constant term, of the form 



5 = b 1 ir 2 9 z + b 2 n 4 9 4 + 0(9 



(5.13) 



From (5.12), we obtain 
1 







A Q - 2h + 



N 



7T 



2 Q 2 



+ A l b l + A Q b\ - 2b 2 



N 



Ac 



1 



12iV 2 12iV 3 



7T 



r 4M 



9 4 + 0(9 b ). (5.14) 
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Setting each term in (5.14) equal to and solving for b\ and b 2 we have the following: 



Proposition 5.1. In the notation above, if < 9 < 1/ir and N is sufficiently large then 
5 = bm 2 9 2 + b 2 n^ + 0(9 6 ) where 



7 A) 1 

h = Y + 2N 



An 



6N 2 24N 3 '' 



(5.15) 
(5.16) 



with Aj given in (5.9). 



Note that in this analysis we have treated F (and hence Aj) as being fixed, in the sense 
that we assume its zeros do not vary with 9. In the next section we will show that when 
f(z) is the characteristic polynomial of a matrix drawn from the CUE, this can be justified 
up to 0{9 7 ). 



Using (5.13) and (5.15), we can determine the distribution of S from the distributions of 



9 and the Aj. For small 9, this comes from the tail of the nearest-neighbor spacing. 

5.2. Nearest-neighbor spacing. Proposition 5.1 provides a formula for 5 = N(l — z'), 



but what we really want is the distribution of 5* := N(l — \z'\). So by (5.4) and (5.13), and 
writing Bj = Re(fy), we have 



5* = B l n 2 9 2 + B 2 ^9^ + I 6» b + — + 



4M 



. 



B 1 tx 2 9 2 + B 2 ^9" + O [ 9 b + 



4M 



N 



N 2 



(5.17) 
(5.18) 



The second line is a corollary of Proposition 8.1, because 5<^9if9<l/n. 



Suppose for the moment that B\ and B 2 were constants (instead of being random). We 
would have 5* = g(9) = B 1 u 2 9 2 + B 2 n 4 9 A + N (9 6 ) where 9 is random with p.d.f. < \3.12} given 
by the nearest neighbor spacing of eigenvalues of unitary matrices. Then the distribution 
function of 5* would be given by 



p 2 {g-\s)) _ B x 



3/2 



6n 
1 

7T 



1 - 



N 2 



,1/2 



-5/2 



45 

fO^(. 5/2 ) 



1 - 



+ 



2N 2 2N A 



5b; 



- 7 ' 2 B 2 



12 



1 - 



N 2 



,3/2 



5.3 



we show that B\ 



(5.19) 
It is 



It turns out that B\ actually is a constant: in Section 
fortunate that Bi is a constant, otherwise it could be difficult to determine the expected 
value of quantities like B l 3 ^ 2 . 

The contribution of B 2 takes a bit more work. If B 2 was independent of 9 then we could 



just average over the possible contributions of B 2 to (5.19). That is, in (5.19) replace B 2 by 



its expected value. This can be computed from the expected values of various combinations 
of Aq and A\. But B 2 is not independent of 9. However, it is independent of 9 to leading 
order. Our specific concern is the distribution of the other roots when 9 is very small. 
This approximates the polynomial having a double zero. Since the next-nearest-neighbor 
spacing of U(iV) eigenvalues vanishes to order 7, the dependence of B 2 on 9 is only to order 
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0{9 7 ). Thus, for the terms we are computing for the distribution of 5* we can treat B 2 as 



independent of 9. The expected value of B 2 is calculated in Section |5.3 

5.3. Expected value of A,. We have 

1 /p'\U) 
A,: ' 



3 ' Ni+ l j\ \F 



N-2 



n=l 

where t\, t 2 , ■ ■ ■ are the arguments of the zeros of F(z). Since 



1 1 % (t\ 

+ — cot I — ) (5.21) 



we see that 



so 



as claimed. Note also that 



1 - e u 2 2 \ 2 

Re(A ) = (5.22) 
B 1 = Re(h) = i (5.23) 

<4>> = \~Jj> ( 5 - 24 ) 

because the imaginary part of the summand is odd. 

For Aj with j > 1, we require a random matrix calculation. The sum in (5.20 ) is dominated 
by the terms where e ltn is close to 1, so one possibility is to determine the level densities of 
the t n . We will find the expected value of Aj by appealing to prior results on averages of 
ratios of characteristic polynomials [3]. 

We assume that f(z) is the characteristic polynomial of a matrix chosen uniformly with 
respect to Haar measure from the unitary group XJ(N). We restrict to those matrices which 
have two eigenvalues very close to 1, and we wish to determine the joint distribution of the 
remaining eigenvalues. This is very similar to the calculations of Duenez [16] and Snaith [T7j 
for the orthogonal group SO (AT). 

First we restrict the measure on the entire ensemble to determine the measure on the 
remaining eigenvalues. Haar measure on U(iV) is given by 

dfi = C \ \ \e itn - e ltm \ 2 dti--- dt N . (5.25) 

l<n<m<N 

Here and following, C is a normalization constant which may vary from line to line, chosen so 
that the measure has total mass 1. Restricting to those matrices which have one eigenvalue 
at 1 is equivalent to rotating (changing variables) to move an eigenvalue to 1. So we can 
also write the measure as 

dm = C Yl \e itn -e ltm \ 2 j J |e i *"-l| 2 dti---dt A r-i. (5.26) 

\<n<m<N~l \<n<N-l 

The set of matrices which have a repeated eigenvalue at 1 has measure zero, so there is no 
canonical way to restrict the measure. However, we are interested in the limiting case of two 
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eigenvalues which are very close together, so we determine the measure by restricting the 



measure (5.26) to have |izv-i| < t, and then let t — > 0. The resulting measure is 

d/i 2 = C j | \e itn -e ltm \ 2 | | \e itn -l\ 4 dt 1 ---dt N _ 2 . (5.27) 



Kn<m<N-2 Kn<N-2 



Let U 2 (N — 2) denote the ensemble of unitary matrices with joint eigenvalue measure fi 2 . 
Then if g = g(e ltl , . . . , e lt "- 2 ) we have 



(g)u 2 (N-2) = CV(2|A(l)| 4 >u(jv- 2) (5.28) 
where the right side is a Haar measure average, A is the characteristic polynomial, and 

Cn = (|A(l)| 4 >u(W (5-29) 



In other words, an expectation involving repeated eigenvalues on U(iV) is equivalent to an 
expectation on U(iV— 2) with an extra factor of the 4 th power of the characteristic polynomial. 
This is the key observation for computing the expected values of the A, because it reduces 
it to the evaluation of known quantities. 
Specifically, let 



A( e -^)A(e-^)A(e-"3)A( e -^)A(e-^)A(e-fe) 
G{a 1 ,a 2 , a 3 , a 4 ; A, /3 2 ; 7i, I2) = A( e -^)A(e-^) ' ^ 



Theorem 4.1 of [3] provides an explicit formula for the expected value of G for A the char- 
acteristic polynomial of Haar distributed matrices on U(A r — 2). The formula is complicated 
so we do not reproduce it here. This is sufficient to determine the expected values of all the 



quantities in (5.15). The calculation requires the assistance of a computer algebra package. 
We now present the answers, which we determined with the help of Mathematica. 
The normalization constant for the measure [i 2 is (the reciprocal of) 



(G(0,0,0,0;0,0;0,0)) u(w _ 2) = (|A(1)| 4 ) U(7V _ 2) 

~ 12 ~ 12" 

= C N \ (5.31) 
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say. As 6 — > we have 



U 2 (N-2) 



-c 



N 



d ai 



«2,7l 



(G(ai,a 2 , 0,0; 0,0; 7i,0)) 



(ai,Q2,7i)=(0,0,0) 
6 



U(JV-2) 



10 10 5 5 



(Ai) 



AV 

A 



C N 1 + 



C/ 2 (AT-2) 



5 



— - -iV+ — 
15 2 15 



ai=0 
11 



(G(ai, 0,0,0; 0,0; 7i,0)) 



U(AT-2) 



7i= Q l 



Thus, 



(A Ai) 



C N [1 



-N 



d 

da>i 



U 2 {N-2) 

d 







7i=ai 



<9«2 



(a 2 ,72)=(0,0) 



(G(ai, a 2 , 0, 0; 0, 0; 71, 7 2 )) U(JV _ 2) 



-iV 2 



13 



30 10 15 

(£? 2 ) = Re(6 2 ) 



-N 



5' 

J_ _ — 
48 ~ 48 



N 



-i 



-2> 



(5.32) 
(5.33) 



Inserting this into (5.19) gives the expansion claimed in Theorem 3.1 

6. The Riemann zeta-function 

We do analogous calculations for derivatives of the Riemann zeta-function and compare 
our results with data. 
We start with 



COO 



i ir'(i s + i) 



i 2r(|s + i) 



1 



1 



s-p p 



(6.1) 



where b = log 2n — 1 — ^7. As in the polynomial case, we assume there are two very closely 
spaced zeros of the (^-function and look for the nearby zero of Suppose the closely spaced 
zeros are 

p ± = l - +l (t±\Q) (6.2) 

with 

(6.3) 



s = - + X + it 
2 



14 DUENEZ, FARMER, FROEHLICH, HUGHES, MEZZADRI, AND PHAN 



a zero of 
Using 



we have 



where 



Note that 



and 



£( a )=log( a ) + 0(l/ a ) (6.4) 

, 4 1, 1 1 v-^ / 1 l\ , . 

= 6 *_ iogt+_ + - +J2 " + " > 6 - 5 

2 s'-p- s'- P+ ^ ± \s ~ P PJ 

b* =b+ 1 -\og(2)-z^ + 0(l/t). (6.6) 

1 1 8X . . 

+ 3 1" = ,^9 , Q 9 ( 6 - 7 ) 



s'-p_ s'-p+ 4X 2 + 6 2 
1 1 = X-i(t-j) |-j7 

V5 , /j. -A9 ^ 1 , .9' V U ' C 



which has real part 



s'-p p X 2 + (t- 7 ) 2 i+7 2 
X 



i 

i 



X 2 + (t- 7 ) 2 + I+ 7 2 ' (6 ' 9) 



As in the polynomial case, we rescale: 

X —x I log t 

e=2rr9/\ogt. (6.10) 
Note that this is analogous to the rescaling in the unitary case because N w log(t/27r). We 

°^" + / 4 logt + J^ + * logi £^7W (6 - n) 

where b** is a real constant, / is purely imaginary, and the sum is over the rescaled zeros 
7* = logt(t — 7)/27r. We follow the same procedure as in the U(iV) case. First multiply 
through (6.11) by (x 2 + 7r 2 # 2 )/logt and expand the final summand as a series in x, giving 

= 2x + (x 2 + ix 2 6 2 ) (-- + xa t - x 3 a 2 + 0(x 5 ) J + smaller terms, (6.12) 



where 

«i = E(SW- (6 ' 13) 

Note that ( |6.12| ) has the same form as Q5.12| ). Now write x = /3i7r 2 2 + /3 2 7r 4 4 + O(0 6 ) and 
gather terms to get 

- ^] n 2 9 2 + ( 2/?2 - \p\ + fact) 7i 2 9 4 + smaller terms. (6.14) 

Thus, 



2 V 2 



h ~ \ (6.15) 
^ 2 ~ 64 ~ l au ( ^ 6 ' 16 ' ) 
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which exactly corresponds to the U(iV) case in Proposition 5.1 

The sum over zeros ot\ is similar to the expression for A\ in the unitary case. We can 
determine the expected value of the sum if we assume CUE statistics for the zeros of the 
Riemann zeta-function, restrict to having two closely spaced zeros, and find the one-level 
density of the remaining zeros. That calculation is in Section[7j In the notation of Lemma 7.1 
we have 



(6.17) 



Note that this is the same as the expected value of A\. 

Thus, assuming that the spacing of zeros of the Riemann zeta function has the same 
distribution as the spacing of eigenvalues of random unitary matrices, we find that the 
leading order behavior of zeros of £' near the |-line is the same as that of zeros of A' near 
the unit circle. 

we show the empirical distribution of the zeros of £' for 10 6 



< t < 10 b 



In Figure 6.1 

60000. The general shape of the distribution shows a striking similarity with the zeros of 
the derivative of characteristic polynomials of unitary matrices. 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 



Figure 6.1. Normalized distribution of the real part of the zeros of C( s )- Data 
is for the approximately 100000 zeros with imaginary part in [10 6 , 10 6 + 60000]. 



7. Calculation of the 1-level density 

We prove the following 

Lemma 7.1. Fix a,b > —1/2. Ifti,t 2 , . . . ,tM ore independently distributed with respect to 
the probability measure 

d/i = = Yl \e itn -e itm \ 2 Yl \e itn -l\ 2a \e itn + l\ 2b dt 1 ---dt M , (7.1) 

M,a,b x< n<m <M \<n<M 

then the large-M limiting (rescaled) 1-level density of the normalized values tj = t.jM/2-K is 
given by 

W[ a ' b \t) = (V,(7rf) 2 + J a+| (vrf) 2 ) - a^J a ^t)J a+l {-nt). (7.2) 
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The measure div- a,b > is Haar measure on U(M + a + b) restricted to those matrices which 
have a eigenvalues equal to 1 and b eigenvalues equal to —1. Note that w[ a ' b ^ is independent 
of b. 

Proof. For fixed a,b > — | define cj(z) — \z — l\ 2a \z + l\ 2b . Let {0 n }^L o be the sequence in 
C[x] uniquely determined by the following requirements: 

(1) For all n > 0, 4> n is of degree n and has positive leading coefficient. 

(2) For all m, n > 0, 

1 



>> ■■= i~ I Me^Me^wie^dt = 5 mn , (7.3) 

where 5 mn is the Kronecker delta. 

Then {<f) n } is the sequence of normalized orthogonal polynomials on the unit circle with 
respect to the measure dv{z) = (2ii)^ 1 uj(z)d£(z) (where di is the arc-length element.) 
Let 

M-l 

K M (z,w) := ^ 4>n{z)<p n {w) (7.4) 

n=0 

be the projection kernel onto polynomials of degree less than M with respect to the inner 
product (7.3). By the Gaudin-Mehta method, the probability measure (7.1), when regarded 
as a measure on the unit circle, can be rewritten as 

1 n 

(note that the normalization constant Cm,u is no longer needed). Then the 1-level measure 
is w[ M \z)dv{z), where 

wi M) (z) = K M (z,z). (7.6) 

(The normalization above is such that the total mass of the 1-level measure is equal to M.) 

Let the "dual" 0* of a polynomial <f>(z) = c n z n + Cn-iz" 1 ' 1 + ■ • • + c\z + Co of degree n be 
the polynomial 

4>*{z) = z^iz- 1 ) = c n + c n _x + ■ ■ ■ + c x z n - x + c z n . (7.7) 
Then we have the following formula of Szego for the projection kernel ([21], Theorem 11.4.2): 



u f \ 0m(^)0m( w ) - 4>m(z)(/>m(w) , . 

K M {Z,W) = z • {l-o) 

1 — zw 

(This formula is analogous to the classical one of Christoffel and Darboux for the projection 
kernel of orthogonal polynomials on the line.) 

In view of (7.6) and (7.8), in order to find the rescaled limit of the 1-level measure as 
M — > oo it will suffice to derive the asymptotic behavior of the orthogonal polynomials (j) n 
as n — > oo near the point z = +1. Theorem 7.2 and formula (7.16) below are the key 
ingredients, but first we need to introduce some notation. 

Denote by P^' b ^ the classical Jacobi polynomials: they are orthogonal in the interval 
[—1,1] with respect to the measure 

w (a ' b \x) := (l-x) a (l + x) b (7.9) 
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and are normalized as follows ([21], Equation 4.3.3): 



"•n 

Let also 



\Pi a ' b \x)\ 2 w {a ' b \x)dx 



2 a+b+i Y(n + a + l)T(n + b + l) 



^(a+1/2,6+1/2) 



2n + a + b + 1 T{n + l)F(n + a + b + 1) ' 



2<i+&t(a-l/2,6-l/2) 



(7.10) 
(7.11) 



Theorem 7.2. EE 



-n+l. 



-1 



2 + 2" 



WiW = cpt 1/2 *- 1/s > 

Letting c n = (a + b)/(n + a + b) , we have 
A 



1\ p(i+l/2,6+l/2) 



n-1 



2 + 2 



-1 



+ D(z-z 



-lx p (a+l/2,6+l/2) / ^ + ^ 1 



n-1 



(7.12) 



'7T 1 + C r . 



2 K 

'7T 1 — Cy, 

2 /i- 



5 



2 
2 



7T 1 — Cr, 



2 /i 



n-1 



' 7T 1 + 



2 ft, 



n-1 



Equation (7.12) appears as (11.5.4) in Szego's book (except for an obvious typograph- 
ical mistake therein.) The constants A,B,C,D can be easily found using Szego's equa- 
tion (11.5.2) together with the fact that 02n-i is a polynomial, which forces the coefficient 
of z~ n on the right-hand side of equation (7.12) to vanish. We omit the details. 



A meaningful rescaling of the 1-level measure is achieved via the change of variables 

2vr£ 



t 



M 



(7.13) 



Indeed, one finds that the limiting 1-level measure (as M — > oo) is Wi(£)d£, where 

WMO = 0, ( 7 - 14 ) 

77) = lim ^K M (e 2 ^ M ,e^ M )Ju(e^/M)uj(e^/My (7 .i 5 ) 

(Actually, all limiting local correlations can be expressed in terms of the limiting kernel K^, 
not just the 1-level density.) 

It remains to compute K 00 (^,r]). It suffices to use equation (7.12) in formula (7.8) and an 
asymptotic formula by Szego's (see [21], equation (8.21.17)): 



(Nt) + t a+ 0(n a ), (7.16) 

valid for fixed a > —1, b £ M, in the range < t < c/n for any fixed constant c > 0, where 
N = n + (a + b + l)/2 and J a is a Bessel function of the first kind. 
A tedious but straightforward computation finally gives: 



Koo(^v) = y^-^f^- (Ja+l/2K)Ja-l/ 2 M " Jar-l/M)J*+l/2(*V)) (7.17) 



*oo(£,0 = \ {<[J a+1/2 (7T0 2 + </a-l/2«) 2 ] " 2aJ a+1/2 «) ^-1/2^) } ■ (7-1* 
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Remark. The kernel K^^^rf) can be used to compute n-level correlations and spacing sta- 
tistics through (matrix or operator) determinants (see [22] for an explanation and proof of 
these applications). Incidentally, for the purposes of evaluating such determinants, the factor 
e «r(»?-0 ma y be suppressed in (7.17) (this is tantamount to conjugating the corresponding 



integral operator by a unitary transformation). 

□ 

8. Appendix: The Domain of Analyticity of the Root of the Derivative 

by Eduardo Duenez 

Define as before e(x) := e 2mx for x real. For any fixed N > 3 define ejv(x) := e(x/N) = 
e 2mx / N . Consider N unit complex numbers 

ejv(0o), e N (-9 ), e N (9 1 ), e N (9 2 ), • • • , e N {9 N - 2 ), (8.1) 

where 

< 9 < 0, < N-6 Q , j = 1,2,..., N -2. (8.2) 
(I. e., the open arc centered at 1 of the unit circle \z\ — 1 having endpoints cn{±9q) contains 



none of the remaining N — 2 numbers). The inequalities (8.2) define a "pyramid" V 



contained in the cube [0, N] N 1 . For fixed < T < N/2, denote by the closed truncation 
of Vn at height T. Then is defined by the inequalities 

9 <9 j <N-9 , j = l,2,...,N-2; 

0<9 <T. (8.3) 

For notational convenience we will set 6 = (9i, 9 2 , ■ ■ ■ , 6*^-2) and Oo = (9 ; 0). Finally, let 

N-2 

f(z) = /(6 ; z) = (z- e N (9 ))(z - e N (-9 )) J] (z - e N (9,)) (8.4) 



be the monic polynomial of degree iV with roots (8.1). 
In this section we prove the following result. 

Proposition 8.1. With the above notation, for every T < - there exists N(T) such that for 

all N > N(T) and for each Bo in V^, the derivative f'(z) = ^/(@o! z) of f has a unique 
root z' = ?"(Oo) in the open disk with diameter [eAr(— #0), e/v(#o)]- Moreover, q is an analytic 

(T) 

function of 0q in the interior ofV N . 



We remark that, by the Gauss-Lucas theorem, the root z' alluded to in Proposition 8.1 
must lie on or to the left of the vertical diameter \e^{— 9 ), e^^o)]. 

Proof. Fix T < 1/tt and consider iV as a parameter (N > 3) for the time being. Let 

N-2 

F(z)=F(e;z):= \{(z-e N (9 3 )), (8.5) 
i=i 

so 

f{z) = {z- e N (9 ))(z - e N (-9 ))F{z) (8.6) 

and 

f -(z)= g ( z ) + L(Q ] z), (8.7) 
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where 



1 



9W 



L(e;z) 



+ 



F' 



-e N (6 ) z 

N-2 

w = £ v. 



1 

- e N ( 
1 



-0o) 



e N {9j)' 



(8.9) 



Define Cn(x) := cos(27rx/A^) and Sjv(x) := sin(27rx/iV), so e^x) = Cn(x) + is^ix). 
Parametrize the boundary of the disk with diameter ejv(=t0o) (i- e., the disk \z — cn(Qq)\ < 
s N (9 Q )) as: 

^(0) = c J v(^o)+^V(^o). (8.10) 

Since roots of / occur at eAr(±#o) whenever / has multiple roots there, it is best to work 
instead with a slightly deformed contour C obtained as the boundary of the "twice bitten" 
disk 



\z - c N (6 
\z - ejv(-0 o )| > e 



< s N (6 ) 



z - e N (9 )\ > e 



missing tiny e-neighborhoods of the points ejv(±0 o )- We still assume that C is parametrized 
by (8.10), except for in a small (^-neighborhood of any multiple of it. (We will not write 
down the exact parametrization of C for such values of since its precise form will not be 
needed.) 

Write C as a union of four pieces: C = Cd U C r U C u U Cg (down, right, up, left), where 





= W) 


n-5<(t)<7T + 5} 


C r 


= {z(4>) 


-71 + 5 < < -5}, 




= {*{<!>) 


-5<<f>< 5}, 




= M4>) 


5<(f)<n-5}. 



The assumed inequalities (8.3) ensure that there are no zeros of / anywhere on C (and, a 
fortiori, no multiple zeros); hence, every zero of /' inside C is a zero of /'// with the same 
multiplicity. It suffices to ensure that that (for all sufficiently small e) the contour C encloses 
exactly one zero of /'//. Note that, since /(0o; z) has no zeros on or inside C, f'/f has no 
poles there either. By the Argument Principle, the claim in Proposition |8.1| regarding the 



uniqueness of the zero z' of /' is equivalent to showing that the image T> of C under /'// 
has index 1 about the origin for large iV and all sufficiently small e. 

We let w{4>) := y(z((j))) and denote the images of the four pieces Cd,C r ,C u ,C£ of C under 
/'// by V d .V r .V u .V e . 

The idea of the proof is very simple. As we shall show, the dominant part of the logarithmic 
derivative f / f(z) (for z in C) is g(z). The image of C under g(z) is easy to describe explicitly; 
indeed, it is a curve £ having index 1 about the origin. The curve T> can be regarded as a 
perturbation of £. We show that, if H and e are sufficiently small, then the remaining terms 
l/(z — ejy(0j)), j = 1, 2, . . . , N — 2, are small enough to keep the index of V equal to that 
off. 



We denote by £ T , £ U} £g the images under g of the respective parts of C (figure 8.1 ). 
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Figure 8.1. The curve £. 



• £ r is parametrized as 

gW)) = 9 , f i„o + 6 < * < -5. (8.11) 

Thus, £ r is a twice-traversed straight segment starting at the faraway point H = 
g(z(--K+5)) = l/(2s(6> )sin<5), moving leftward to the point A = g(z(0)) = l/(2s(0 o )) 
and retracing itself back to B = g(z(—5)) (here B = H). 

• £ u is a large arc BCD on the upper half-plane. (£ u is essentially a semicircle, for e 
small. This is easily seen from the fact that the dominant term in f'/f(z) for z very 
close to eAr(# ) is l/(z — e N(9 )), namely an inversion with center ejv(#o)> taking the 
tiny (almost) semicircle C u to a huge (almost) semicircle £ u .) £ u escapes any bounded 
region of the plane as e approaches zero. 

• Ei and Sd are obtained from £ T and £ u through central symmetry with respect to the 
origin. 

It is clear that £ has index 1 about the origin for all sufficiently small e. 

For notational convenience, denote by A = w(—ir/2), B = w(—5), . . . , H = w(—tt + 5) the 
points on V analogous to A, B, . . . , H on £. 

(T) 

We claim that if N is large enough and 6o € V N , then the index of V about the origin 
is also 1 for all sufficiently small e and all iV > N(T). Note that e is allowed to depend on 
©o, whereas the lower bound N(T) for iV depends only on T. 

For notational simplicity, we set 0jv-i := — #o- Write 

L( \ = _i_ \^ 1 

/ W ^-e w (^ ) + ^z- ew (^)' 

i2(e ;^) 

where m > 1 is the multiplicity of the zero ex^Oo) of /(z) and 9j, j — m, . . . , N — 1, are the 
remaining zeros of f(z). 

Let e > be the minimum of the distances from eAr(6' ) to the other ejv(^), m < 
j < N - 1. Now let e = e /iV. If \z - e N (9 )\ = e, then \z - e N (9j)\ > e (m < j < 
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iV — 1), so |i2(0o;z)| < (N — m)/eo- Moreover, m/(z — ejv(^o)) maps C u into (part of) 
the upper half-circle S: 



\w\ 



on the interval 
our claim that T> n only intersects the 



m/e, Hew > 0. Therefore, D u will be supported on the 
^-^-neighborhood T of S. Since T intersects the imaginary axis Hew 

■ \m _ N-m m , N-m ] _ ■ \ m+N(m-l) N+m(N-l) 

positive imaginary axis is proved. The claim that T>i only intersects the negative imaginary 
axis for suitably small e is proved along identical lines. 

In order to prove the claim that D has index 1 about the origin, it suffices now to show 
that T> T is contained in the right half-plane Re(u>) > and T>i in the left half-plane Re(w) < 
0. The truth of said statement for X> r is obvious, since each of the terms l/{z — e^iOj)) 
(J = 0, . . . , N — 1) in f'/f(z) (cf., equations (8.7)-( pT9 )) has positive real part for z on C r . 



Now note that 



5<6<tt-8, 



(8.12) 



2sn{0o) sin< 

is real and positive; we anticipate it to be the dominant term of the logarithmic derivative 
— /'/ f(z). We will show that the real parts of each of the terms z(6 )_l,^( .) (for 5 < <fi < n — 5 
and j 



1, . . . , N — 2) are bounded above by a suitably small multiple of —g(z(<f>)). 



Lemma 8.2. For all z with \z\ < 1 we have 

max Re 

ICI=i z - C 



Re(z) 



1 - \z\ 2 

This lemma can be proved by rewriting the equation |£| 

1 



(8.13) 



1 in the form 



w + 



1 



in terms of the variable 



w 



whence the result follows trivially. 

Lemma 8.3. For N > 3, < 6q < I, < <p < it and all xj) we have 



Re 



z(4>) - e N (ip) 



< 



n9o 
N 



2^(7 + 4^)^)^(0)), 



(8.14) 



First, since —g(z ((/))) > for < <\> < it (cf., equation (8.12)), it suffices to prove the 
upper bound 

2^(7 + 4^)11 



TWo 

N 



for the quantity 

h(9 ;<f>,ip) 



Re 



By Lemma 8.2 



g{z{4>)) 

h{6 ]<f>,il>) < 



- Re(z(0)) 



AT2 



sn(@o) sin cj) Re 



z((j)) - e N (i/j) 



115) 



(8.16) 



s N (V ) sin. 
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From Equation (8.10) it quickly follows that 

2 = sn(29q) sin0, and 



Re(>(0)) = c N (9 Q ) - s N (9 ) sin0. 



Thus, 



h(9 -A,tp) < 



1 - c N (9 ) + s N (9 ) sin < 



tan 



s N (29 ) sin 
2n9 



< - tan 
~ 2 



JV 
2ti9 

~aT 



smra + 



sec 



sec 

2tt#o 
~N~ 



s N {U )sw.( 
2<k9q 



1 - c N (9 ) + s N (9 ) sin, 
2c N (9 ) 



(8.17) 



Let us now use the Taylor formula? with remainder 

| tan 9 — 9 | < — sup | tan" d \ 

2 \q\<o 

| sec 9 — 1| < — sup | sec" i?| 
2 |^|< 

with := 2tt9o/N. Both | tan"$| and | sec"??! are even functions of increasing with |$|, so 
their respective suprema are bounded by | tan"(27r/3)| = 14 and | sec"(2vr/3)| = 8^ since 
•d = 9 = 2n9 /N < 2tt/3, by the assumptions 9o < 1 and JV > 3). The inequalities 

tan9 < 9 + 79 2 , and 
sec9 - 1 < Ay/ZO 2 , 



follow immediately. These inequalities together with (8.17) complete the proof of Lemma 8.3 
upon setting 9 = 2tt9 /N. 



To complete the proof that T>i is contained in Re w < 0, it remains to show that 

7', 



- Re 



/ 



for S < 6 < n — 5. 



(8.18) 



But 
-Re 



*(</>)) 



JV-2 

-^(0))-^Re 



> -tfW)) 
= -<K*(0) 



1- (JV-2) 



z(4>) - e A r(6' j 

7r(9 



since —g(z((f>)) > 



by Lemma 8.3 



a a% eft 



say, for suitable absolute constants U, V, W. As JV — > oo, the last bracketed quantity above 
has the limit 1 — 7i9q. As long as \9q\ < T < there will exist JV(T) such that said quantity 
is positive for JV > JV(T). This completes the proof of (8.18) and of the uniqueness of the 
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desired root z' = <r(@o) of /'. The analyticity of ? as a function of Bo follows from the joint 
analyticity of the function /(Bo; z) in the variables 0q and z together with the formula 



1 

2vri 



*/"(e ; z) 



dz. 



119) 



/'(e„;z) 

While C depends on the choice of a fixed e, it is clear that e can be chosen so that z' = ?(0o) 
is uniquely defined by the integral (8.19) for all O in any desired compact subset of the 

(T) (T) 

interior of V N . This is enough to ensure that q is analytic in the whole interior of V N and 
concludes the proof of Proposition 8.1. □ 
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